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Q \ Abstract 

We observed outer regions of a bright cluster of galaxies A2142 with Suzaku. Temperature and brightness 
_J ■ structures were measured out to the virial radius (r2oo) with good sensitivity. We confirmed the temperature drop 

from 9 keV around the cluster center to about 3.5 keV at r2oo, with the density profile well approximated by the /3 
model with j3 = 0.85. Within 0.4 r^oo, the entropy profile agrees with r • , as predicted by the accretion shock model. 
The entropy slope becomes flatter in the outer region and negative around r^m. These features suggest that the 
> , intracluster medium in the outer region is out of thermal equilibrium. Since the relaxation timescale of electron-ion 

Coulomb collision is expected to be longer than the elapsed time after shock heating at r2oo, one plausible reason 
of the low entropy is the low electron temperature compared to that of ions. Other possible explanations would 
be gas dumpiness, turbulence and bulk motions of ICM. We also searched for a warm-hot intergalactic medium 
around r^ao and set an upper limit on the oxygen line intensity. Assuming a line-of-sight depth of 2 Mpc and oxygen 
abundance of 0. 1 solar, the upper limit of an overdensity is calculated to be 280 or 380, depending on the foreground 
ly-j \ assumption. 

\Q . Keywords: galaxies: clusters: individual (Abell 2142) — X-rays: galaxies: clusters — X-rays: ICM, WHIM 

lO ' (Warm Hot Intergaractic medium) 
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1. Introduction cluster outskirts are connected to the surrounding large-scale 

structure, where the gas is falling towards the cluster potential 
The Cold Dark Matter scenario of the cosmic structure and possibly subject to shock heating. Once disturbed by the 
formation predicts that clusters of galaxies are formed via subcluster collisions, it should take long time for the gas to set- 
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collisions and mergers of smaller groups and clusters. As tie due to the low density and large spatial size. Therefore, the 

shown by numerical simulations, merging plays a critical role outskirts of merging clusters offer us opportunities to look into 

in the cluster evolution. X-ray observations have provided the gas in its transition to the thermal and ionization equilibri- 

many pieces of evidences for cluster mergers, through imag- ums. 

ing infalling subclusters and disturbed, irregular morphologies Recent X-ray studies with Suzaku showed temperature 
of intracluster medium (ICM) (e.g., Forman & Jones 1982). structure of ICM to the virial radius (r2oo) for several relaxed 
Furthermore, even in apparently relaxed clusters, discontin- clusters (George et al. 2008; Reiprich et al. 2009; Bautz et al. 
uous ICM structures so-called "cold fronts" are sometimes 2009; Hoshino et al. 2010; Simionescu et al. 2011). The re- 
found from the Chandra observations. The cold fronts are in- suits show a systematic drop of temperature by a factor of ~ 3 
terpreted as contact discontinuity caused by the subcluster col- from the center to r^oo- Also, the entropy shows a flattening or a 
lisions (Markevitch & Vikhlinin 2007). Although the morpho- small drop, after a monotonous increase with radius, around the 
logical data are accumulating, the physical states of the ICM, outermost region. It has been interpreted as the deviation be- 
particularly the ionization equilibrium state and the kinematics, tween electron and ion temperatures in this region. Numerical 
are yet to be clarified. simulations for relaxed clusters also suggest that the outskirts 
In the study of the dynamical evolution of clusters, we will of clusters are not in hydrostatic equilibrium (e.g. Burns et al. 
focus on the cluster outer regions since they are expected to 2010). 

contain important information on the formation process. The Compared with the relaxed clusters, however, there is lit- 
tle information about the ICM properties in the outer regions 



Last update: 2011 June 30 c , , T , • n .u 

v of merging clusters. In those dynamically-young systems, the 
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ICM are thought to be in the early stage of thermal relaxation. 
Thus, the X-ray observations of merging clusters and compar- 
isons with the relaxed clusters will give us new, valuable infor- 
mation on the ICM properties at large radii and help us to draw 
the global picture of formation and evolution of clusters. 

The surface brightness around the cluster virial radius is 
much lower than the central region (typically, by a factor of 
10 4 ). Thus, the detailed estimation of foreground and back- 
ground emission as well as careful assessment of all the sys- 
tematic errors is important in this kind of study. The XIS in- 
strument on Suzaku has a low and stable background and the 
behavior of the non-X-ray background (NXB) is well known 
(Tawa et al. 2008; Koyama et al. 2007), which makes XIS the 
most suitable for the study of the cluster outer regions. 

An additional science at cluster outer regions is the search 
for the warm-hot intergalactic medium (WHIM). The filamen- 
tary structure of the local universe has been probed mainly 
through galaxy distribution (e.g. Eisenstein et al. 2005). Its 
detailed structure would be directly observed by the WHIM, 
which has temperature of 10 6 - 10 7 K and contains more than 
half of baryons in the local universe (e.g. Cen & Ostriker 2006). 
Because of its extreme faintness, detailed observations of the 
WHIM will be the subject of future high-resolution X-ray stud- 
ies. The gas density and chemical composition of the WHIM 
are still poorly known both theoretically and observationally. 
The cluster outskirts are the regions where the ICM is con- 
nected to the WHIM and will enable us to place observational 
constraints about the WHIM properties. As shown, by e.g., 
Takei et al. (2008), Suzaku XIS is able to set strong constraints 
about the redshifted oxygen emission. 

Abell 2142 (A2142, z = 0.0909) is a bright cluster of galax- 
ies, having a high ICM temperature of kT » 9 keV. This ob- 
ject is also known as the first cluster in which the cold fronts 
have been detected (Markevitch et al. 2000). The cold fronts in 
the south and the northwest are 0.'7 (or 70 kpc) and 2.'7 (or 
270 kpc) off of the cluster center, respectively, and a sharp 
surface-brightness drop by a factor of about 2 is seen. Since the 
temperature and density distributions suggest that the pressure 
is constant across the cold front, it is considered as a contact 
discontinuity. The presence of these structures naturally indi- 
cates that A2142 is a merger and the subcluster infall seems to 
be occurring along the northwest-southeast direction. Actually 
the overall cluster emission is elongated in this direction. This 
is likely to coincide with the large-scale structure, and the mat- 
ter density seems to be enhanced along the filament. 

A2142 is also suitable for the search of emission lines from 
the WHIM. Given the cluster redshift of ~ 0.1, the OVII line 
(the rest-frame energy of 0.65 keV) is shifted to 0.57 keV and 
then falls into a gap between the Galactic OVII line energy and 
the instrumental line features. Therefore, it makes it possible 
to distinguish the WHIM emission from the local one. In addi- 
tion, because the redshift of A2142 is not too high, the oxygen 
line from the WHIM can be measured with good sensitivity. 
With these purposes, we have carried out Suzaku observations 
of the northwest offset regions of A2142. 

In this paper, we will use Ho = 70 km s - ' Mpc -1 , Q M 
0.27 and £2a = 0.73. This cosmology leads to 100.4 kpc per 
arcmin at z = 0.0909. The virial radius approximated by /-200 = 
2.77/!7 ( }«r>/10keV) 1/2 Mpc/,E(z), with E(z) = (fi M (l +z) 3 + 1 - 



Dm) 1 ' 2 (Henry et al. 2009). For our cosmology and redshift, 
r 20 o is 2.48 Mpc (= 24 .'8) with kT = 8.7 keV. In this paper, 
we adopted the solar abundance defined by Anders & Grevesse 
(1989) and the Galactic N H by Dickey & Lockman (1990). The 
errors are in the 90% confidence for a single parameter. 

2. Observations and Data Reduction 

As shown in figure 1, we performed four pointing observa- 
tions in 2007 with the XIS instrument around Abell 2142. The 
central pointing observation was performed in June, and the 
other three in September to October. The cold front feature 
indicates that there is an ongoing merger in the north-west to 
south-east direction. This suggests that matter would be falling 
in along this merger axis, possibly from a large-scale filament. 
We allocated the observed regions to be successively offset to- 
wards the north-west direction. The regions are designated as 
Center, Offsetl, Offset2, and Offset3 with the exposure times 
51.4 ks, 37 ks, 58 ks, and 24 ks, respectively. The observation 
log is shown in table 1. The outermost observation reaches 
twice the virial radius (49.'2 ~ 4.92 Mpc) from the cluster cen- 
ter, in which we planned to search for emission from the warm- 
hot intergalactic medium (WHIM). 

Three out of the four CCD chips were available in these 
observations: XIS0, XIS1 and XIS3. The XIS1 is a back- 
illuminated chip with high sensitivity in the soft X-ray energy 
range. The effective area has been siginificantly reduced by the 
time of the observations due to the contamination building up 
on the IR/UV blocking filters. This effect along with its un- 
certainty are included in the effective area in our analysis. We 
used HEAsoft ver 6.9 and CALDB 2010-12-06. All the XIS 
sensors were in the normal clocking mode, and the spaced-row 
charge injection (SCI) was applied. 

We extracted pulse-height spectra in 9 annular regions with 
boundaries at 0', 2.'5, 5.'0, 7.'8, 10.'3, 12.'9, 16.'8, 24.'6, 31.'1, 
38.'8 centered on (15 h 58 m 16. s 13, 27°13'28") from all the XIS 
events. We analyzed the spectra in the 0.5-10 keV range for 
the FI detectors and 0.35-8 keV for the BI detector. In all an- 
nuli, the calibration source positions were masked out using the 
calmask calibration database (CALDB) file. 

3. Background Analysis 

In the present study of the ICM emission in the cluster 
outskirts, correct estimation of the background is of utmost 
importance. As a standard practice, we assume three back- 
ground components: non-X-ray background (NXB), cosmic 
X-ray background (CXB) and Galactic emission. The Galactic 
emission component has a spatially variable spectrum. We es- 
timated the Galactic background spectra using two Suzaku ob- 
servation data, one was Offset3 and another one was observa- 
tion of TCrB (Observation ID = 401043010) which was located 
at 1 degree south of A2142. We will show each background 
component in this section. 

3.1. Cosmic X-ray Background 

Since the CXB consists of many extragalactic point sources, 
we tried to remove the sources as many as possible and 
then modeled the remaining emission with a power-law. We 
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Table 1. Log of Suzaku observations of Abell 2142 
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Position(Obs. ID) 


Start 




End 


Exp. time (ks)* 


Exp. time (ks) 1 


Center (801055010) 


2007 


Jun 


04 


2007 


Jun 


05 


51.4 


45.2 


Offsetl (802030010) 


2007 


Oct 


04 


2007 


Oct 


05 


37.6 


26.6 


Offset2 (802031010) 


2007 


Sep 


15 


2007 


Sep 


17 


57.6 


41.3 


Offset3 (802032010) 


2007 


Oct 


29 


2007 


Oct 


30 


23.7 


20.4 


TCrB (401043010) 


2006 


Sep 


06 


2006 


Sep 


08 


46.3 


36.8 



*:COR2 > GV t COR2 > 8 GV 




Fig. 1. (a):Rosat All Sky Survey image around A2142. Cyan circles show XMM-Newton FOV, and Magenta boxes show Suzaku Observations. White 
contours show galaxy distribution associated with the A2142 taken SDSS catalogue. (b):NXB subtracted Suzaku FI+BI image of A2142 in 0.5-8.0 keV 
band smoothed by a 2-dimensional gaussian with tr = 16 pixel =17" . The image is corrected for exposure time but not for vignetting. Large white circles 
indicate the regions used for spectrum analysis. Large green circle show the virial radius of A2142 (~ 2.5 Mpc). Small white circles show the excluded 
point sources. 



estimated the CXB surface brightness to be 5.97 xlO 
erg cirT 2 s _1 sr~' after the source subtraction, based on ASCA 
GIS measurements (Kushino et al. 2002). We carefully sub- 
tracted point sources brighter than 8 x 10~ 14 erg crrT 2 s , while 
the flux limit of the point sources eliminated in Kushino et al. 
(2002) is 2 x 10" 13 erg crrT 2 s _1 . Our flux limit is thus suffi- 
ciently lower than Kushino et al. (2002). The details of point- 
source subtraction are described in Appendix 1 . 

To estimate the amplitude of the CXB fluctuation, we scaled 
the fluctuation measured with Ginga (Hayashida et al. 1989) 
to our flux limit and field of view, following the analysis by 
Hoshino et al. (2010). The fluctuation amplitude scales as 
(^e.Suzaku/^e.Ginga) -0 ' 5 , with Q e ,Suzaku and Q e>G in g a the effective 
field of views (FOVs) of Suzaku and Ginga instruments, re- 
spectively. We show the resultant relative fluctuation cr// CXB 
for each annular region in table 3, where <x is the standard de- 
viation of the CXB intensity /cxb- 

3.2. Non X-ray Background 

The non X-ray background (NXB) spectra were estimated 
from the database of Suzaku night-earth observations using the 
procedure of Tawa et al. (2008). We accumulated the data 
for the same detector area and the same distribution of COR2 
as the A2142 observations, using an FTOOL xisnxbgen. The 



night-earth data cover 150 days before and after the period 
of A2142 observations. To increase the signal-to-noise ratio 
by keeping the NXB count rate low, we selected durations in 
which COR2 is > 8 GV. The systematic error due to the NXB 
uncertainty was estimated by varying the NXB intensity by 
+3% as in Tawa et al. (2008). 

3.3. Galactic Components 

To estimate the Galactic emission, we examined the spectra 
from two Suzaku observations: Offset3, which showed neg- 
ligible ICM contribution, and TCrB at 1° south of A2142. 
We employed ancillary response files (ARFs) for a spatially 
uniform source filling the FOV. A power-law is used to 
model the CXB in both spectra. The Galactic emission is 
represented as a two-temperature model consisting of an un- 
absorbed ~0.1 keV plasma (LHB; representing the local hot 
bubble and the solar wind charge exchange) and an absorbed 
~0.3 keV plasma (MWH; representing the Milky Way halo): 
apec x + wabs x (apec 2 + powerlaw). The redshift and abun- 
dance of both the apec components were fixed at and unity, 
respectively. 

The best-fit parameters for the Offset3 and TCrB regions are 
summarized in table 2. The temperatures of the LHB and the 
MWH are 0.09 ± 0.02 keV and 0.28 ± 0.04 keV, respectively, 
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Fig. 2. NXB-subtracted XIS BI spectra (Black cross) of Offset3 (left) and TCrB(right), plotted with estimated CXB (black line) and NXB (gray cross) 
spectrum described in Sec 3.1, Sec 3.2 respectively. 



for Offset3, and 0.09 + 0.02 keV and 0.29 + 0.03 keV, respec- 
tively, for TCrB background. The temperatures and intensities 
are consistent with the typical Galactic emission. The observed 
fluxes in the two regions differ by 8%, higher in Offset3. We 
assume that the difference of the two fluxes represents the typ- 
ical uncertainty range and the systematic error of the Galactic 
emission. 

3.4. Background Fraction in Each Region 

Table 3 summarize the information in each annular region 
we analyzed. The columns indicate: the annular boundaries; 
O e , the solid angle of observed areas; Coverage, the cover- 
age fraction of each annulus, which is the ratio of £2 e to the 
total solid angle of the annulus; SOURCEJtATIOJIEG, the 
fraction of the simulated cluster photons that fall in the region 
compared with the total photons generated in the entire simu- 
lated cluster; <t// CX b, the CXB fluctuation due to unresolved 
point sources; OBS, the observed counts; the estimated counts 
for the three background components, i.e., NXB, CXB, and 
the Galactic emission; and the fraction of background photons 
given by / B gd = (NXB+CXB+Galactic)/OBS. 

The NXB counts are calculated from the night earth data. We 
simulated the spectra of the Galactic and CXB components, 
using xissim (Ishisaki et al. 2007) with the flux and spectral 
parameters given in the row of "A2142 OFFSET3 nominal" in 
table 2, assuming a uniform surface brightness that fills the XIS 
FOV. We plot the NXB and CXB BI spectra compared with the 
observed spectra in the background regions in figure 2. These 
simulated spectra gave the counts shown in table 3. 

4. Stray Light 

Stray light is photons entering from outside of the FOV. 
The Suzaku optics often show non-negligible stray light from 
nearby bright X-ray sources (Serlemitsos et al. 2007). The ex- 
tended telescope point spread function also makes contamina- 
tion of photons from a nearby sky. We estimated the contami- 
nation flux based on a ray tracing simulation xissim. The result 
of our simulation is shown in table 4. For each observed an- 
nulus, we calculated the fraction originated from each sky area 



(i.e., annulus). We can see that the "on-source" fractions are 
naturally higher for the central pointing than for the other re- 
gions, since the on-source flux itself is higher. At the same 
time, 20-30% of the detected photons are due to the contami- 
nation, and they are mostly from adjacent regions. Therefore, 
the results indicate that even though the flux contamination is 
not negligible, the origin is limited to nearby regions. We pro- 
ceeded to the spectral analysis without making correction for 
the stray light. 

5. Spectral Analysis 

5.1. Spatial and Spectral Responses 

We need to calculate the spatial and spectral responses for 
the analysis of A2142 data. The response functions for ex- 
tended sources are complicated because they depend on the 
surface brightness distribution of the source. They need to be 
calculated for each annular region. Monte Carlo simulator xis- 
sim incorporates the responses of the X-ray telescope and XIS 
instrument. The ARF generator using this simulator is called 
xissimarfgen (Ishisaki et al. 2007). We used version 2008-04- 
05 of the simulator. The surface brightness distribution is one 
of the input parameters necessary to run xissim and xissimarf- 
gen. Because of the extended PSF, the local efficiency is related 
with the relative flux among adjacent spatial elements. We used 
the /3-model (B = 0.85, r c = 4.'5) based on the ROSAT PSPC re- 
sult as the input X-ray image (Henry & Briel 1996). 

We created ARFs assuming that the input image does not 
vary with energy. The effect of the contamination on the XIS 
IR/UV blocking filter is included in the ARFs based on the 
calibration in November 2006. The normalization of the ARF 
is defined such that the flux given as a result of the spectral fit 
is equal to the entire flux for a given spatial region. 

5.2. Spectral Fit 

We carried out spectral fitting to the pulse-height data of 
each annular region separately. The NXB component was sub- 
tracted before the fit, and the fitting model included the LHB, 
MWH, CXB and ICM components. The spectra from the BI 
and FI sensors were jointly fitted with the same model by min- 



No. ] 



X-Ray Study of the Outer Region of Abell 2142 with Suzaku 

Table 2. Best-fit values of background spectra fitting 



A2142 OFFSET3 



l[0.4-10.0keVl 



' 2|0.4-10.0keVl 



Unabsorbed (keV) 

— o 090 +O020 — 

-MIS 
090 +0 " 13 

°- 095 QQ21 
092 + " : ° 22 

094 +5 ' 5T3 



norm. 



4.90 +4/b 
6 13+^ 

t 97+2.23 
-1-67 
4 01 +6 ' 5 

7.01 



Absorbed (keV) norm! 



n 97«+u.u4i 

279 + °°^ 
-P.-Q33 
281 +0 -°*o 

0.266+M 



51 +0 - 2u 

55 +il22 
Q 12 

q_ 4 7+0.1§ 

53 + " : ^ 
-Q-J3 

)LzL=sm 



nominal 

CONTAMI+10% 
CONTAMI-10% 
NXB+3%+MAXCXB 
NXB-3%+MINCXB 



4.17; 



0.59+0.02 
0.68+0.02 
0.53+0.02 
0.63+0.02 
0.59+0.02 



0.286 



-8:811 

-0.041 



0.68+0.02 
0.74+0.02 
0.63+0.02 
0.72+0.02 
0.66+0.02 



TCrB 



Unabsorbed plasma 



Absorbed plasma 





kT (keV) 


norm\ 


<? 1 " 

° 1T0.4-I0.0ke VI 


kT (keV) 


norm* 2 


9* 

J 2r0.4-10.0ke VI 


nominal 


0.089™ 


3.85«-it 
4.48 + °f° 
3 36 + ^ 2 
4.24 + ^I 

— 1.00 
3 95+5.43 


0.41+0.01 


9Q9 +0.032 

-Q-P20 


52 +u tw 

-Q1Q 


0.72+0.02 


COMTAMI+10% 


0.088!™ 


0.35+0.01 


0291 +H.83¥ 

-Q-P22 


51+0.69 
-Q-Q2 


0.71+0.02 


CONTAMI-10% 

NXB+3%+MAXCXB 

NXB-3%+MINCXB 


0-087-II 
0.086! !j 


0.42+0.01 
0.44+0.01 
0.34+0.01 


n 9OC+0.050 

O.285_ 0Q23 

270 +0 ° 42 
293+0.030 


045 +0.0? 

U - HJ -0.12 

53 +01 ° 

—0.13 

52 +009 

U - J -0.09 


0.72+0.02 
0.70+0.02 
0.60+0.02 



*: Normalization of the apec component scaled with a factor 1/400/r assumed in the uniform-sky ARF calculation (circle radius r=20'). 

Norm= OTr/ «eiHdV/(4ff(l + z 2 )D^) x 10~ 20 cm" 5 arcmin -2 , where Da is the angular diameter distance to the source, 
t: 10~ 6 photons cm _2 s _1 arcmin -2 . Energy band is 0.4 -10.0 keV. 



Table 3. Estimation of CXB fluctuation and Backgrand Count 



Region 




Coverage^ 


SOURCE. 








FI Count (0.5-10 keV) 






(arcmin 2 ) 


(%) 


Ratiojleg^ 


(%) 


OBS (xlO 2 ) 


NXB (XlO 2 ) 


CXB (XlO 2 ) 


Galactic (xlO 2 ) 


fBGDi.%) 


0'-2.'5 


21.0 


100.0 


27.0 


24.1 


1701.8 ±4.1 


4.1 ±0.1 


3.9 ± 0.9 


1.5 ±0.1 


0.6 ±0.1 


2.'5-5.'0 


58.8 


95.4 


31.1 


14.4 


1131.2 ± 3.4 


11.5 ±0.3 


9.9 ± 1.4 


4.4 ± 0.2 


2.3 ± 0.4 


5.'0-7.'8 


96.3 


93.6 


16.0 


11.3 


326.8 ± 1.8 


10.8 ± 0.3 


9.8 ± 1.1 


3.8 ± 0.2 


7.5 ± 1.4 


7.'8-10.'3 


126.6 


86.2 


7.7 


9.8 


122.0 ± 1.1 


9.2 ± 0.3 


7.2 ± 0.7 


3.1 ±0.2 


16.1 ±3.1 


10.'3-12.'9 


128.9 


67.6 


3.4 


9.7 


36.9 ± 0.8 


6.1 ±0.2 


5.7 ± 0.6 


2.5 ± 0.2 


38.7 ±8.2 


12.'9-16.'8 


101.8 


42.5 


1.1 


10.3 


31.8 ±0.6 


7.8 ±0.2 


7.1 ±0.8 


2.8 ± 0.2 


55.7 ± 11.2 


16.' 8 -24.' 6 


242.9 


22.2 


0.8 


7.2 


21.7 ±0.5 


8.5 ±0.3 


5.8 ± 0.4 


2.4 ± 0.2 


77.2 ± 15.6 


24.'6-31.'l 


191.6 


15.8 


0.2 


8.1 


38.9 ± 0.6 


17.6 ± 0.5 


16.8 ± 1.4 


6.4 ± 0.2 


104.9 ± 17.1 


31.'l-38.'8 


241.7 


14.1 


0.1 


7.2 


31.2 ±0.6 


16.5 ± 0.5 


11.9 ±0.9 


5.2 ± 0.2 


107.8 ± 18.4 


Region 


or. 


Coverage"* 


S ource. 








BI Count (0.5-8 


keV) 






(arcmin 2 ) 


(%) 


Ratiojteg^ 


(%) 


OBS (XlO 2 ) 


NXB (XlO 2 ) 


CXB (XlO 2 ) 


Galactic (xlO 2 ) 


f B GD(%) 


0'-2.'5 


21.0 


100.0 


27.0 


24.1 


1042.6 ± 3.2 


3.0 ±0.1 


1.9 ±0.5 


1.2 ±0.1 


0.6 ± 0.2 


2.'5-5.'0 


58.8 


95.4 


31.1 


14.4 


732.5 ± 2.7 


8.2 ± 0.2 


5.1 ±0.7 


3.8 ± 0.2 


2.3 ± 0.5 


5.'0-7.'8 


95.8 


92.5 


15.5 


11.5 


245.4 ± 1.6 


9.8 ± 0.3 


6.1 ±0.7 


4.6 ± 0.2 


8.4 ± 1.6 


7.'8-10.'3 


122.6 


86.2 


7.7 


9.8 


88.8 ± 0.9 


7.8 ±0.2 


4.9 ± 0.5 


3.2 ± 0.2 


17.8 ±3.7 


10.'3-12.'9 


125.9 


64.3 


3.0 


9.4 


30.5 ± 0.7 


4.8 ±0.1 


3.3 ± 0.3 


2.2 ± 0.2 


33.8 ±7.8 


12.'9-16.'8 


105.8 


46.7 


1.2 


10.8 


26.2 ± 0.5 


6.0 ± 0.2 


4.1 ±0.5 


2.8 ± 0.2 


49.2 ± 10.7 


16.' 8 -24.' 6 


233.9 


20.8 


0.7 


8.2 


18.6 ± 0.4 


7.2 ± 0.2 


3.6 ± 0.3 


2.7 ± 0.2 


72.8 ± 15.6 


24.'6-31.'l 


181.6 


13.8 


0.2 


9.3 


31.6 ±0.6 


14.6 ± 0.4 


9.2 ± 0.7 


6.2 ± 0.2 


94.9 ± 16.8 


31.'l-38.'8 


225.7 


11.1 


0.1 


7.6 


24.6 ± 0.5 


12.7 ± 0.4 


6.8 ± 0.5 


5.0 ± 0.2 


99.3 ± 18.4 



*: Solid angle of each observed region, 
f: Fraction of each area to entire annulus. 

t: Fraction of the simulated cluster photons which fall in the region compared with the total photons generated in the entire simulated cluster. 
§: CXB fluctuation due to unresolved point sources. 



imizing the total x 2 value. To increase the signal to noise ra- 
tio of A2142, we used energy ranges of 0.35-8 keV for BI 
and 0.5-10 keV for FT. The relative normalization between 
the two sensors was a free parameter in this fit to compen- 
sate for the cross-calibration errors. The photon index and 
normalization of the CXB, the temperatures and normaliza- 
tion of the LHB and MWH were fixed at the values in ta- 
ble 2. Metal abundances of the LHB and MWH components 
were set to be unity. The Galactic absorption column density 



was fixed at /V H = 4.2 x 10 20 crrT 2 (Dickey & Lockman 1990). 
We checked the influence of an uncertainty in A^h by using an- 
other column density from Leiden/ Argentine/Bonn (LAB) sur- 
vey (Nh = 3.8 x 10 20 cm" 2 ; Kalberla et al. 2005) to confirm that 
the two spectrum fits do not show significant difference. The 
fits were carried out with XSPEC verl2.4.0ao. In the central 
regions, free parameters were the temperature, normalization, 
metal abundance of the ICM component. In the outer regions, 
we fixed the metal abundance of the ICM at 0.2, which is the 
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Table 4. Stray light contamination for the central, OFFSET 1 and OFFSET2 pointings 



[Vol. 



CENTER 



Emission Weighted 


Detector 


0'-2.'5 


2.'5-5.'0 


5.'0-7.'8 


7.'8-10.'3 


10.'3-12.'9 


Radius (arcmin) 


/Sky 












1 


0'-2.'5 


81.5 % 


17.9 % 


0.5 % 


0.1 % 


0.0% 


4^ 


2.'5-5.'0 


26.5 % 


66.9 % 


6.5 % 


0.2% 


0.0% 


59=** 
84^ 


5.'0-7.'8 


5.3 % 


33.8 % 


56.0 % 


4.8% 


0.1 % 


7.'8-10.'3 


4.3 % 


8.1 % 


31.4% 


51.9 % 


4.3 % 


110^ 


10.'3-12.'9 


2.9 % 


6.0% 


6.8 % 


29.1 % 


55.2 % 


OFFSET 1 and OFFSET2 


Emission Weighted 


Detector 


<12.'9 


12.'9-16.'8 


16.'8-24.'6 


24.'6-31.'l 


31.'l-38.'8 


Radius (arcmin) 


/Sky 














<12.'9 


94.1% 


4.7 % 


0.9% 


0.2% 


0.0% 


U.0 +2S 
18.6 +: 


12.'9-16.'8 
16.'8-24.'6 


21.2% 
1.8% 


61.7% 
17.8% 


15.8% 
67.4% 


1.2% 
13.0% 


0.0% 
0.0% 


26 - 3 18 


24:6-3 l.'l 


1.3% 


2.4% 


32.6% 


63.1% 


0.5% 




31.'l-38.'8 


6.7% 


0.0% 


20.0% 


13.3% 


60.0% 



lowest value observed in the outskirts of clusters (Fujita et al. 
2008). 

Figure 3 shows the results of the spectral fit for all the annu- 
lar regions. The parameters and x 1 values are listed in tables 
5 and 6. We obtained fairly good fits for all the regions with 
reduced^ 2 values less than 1.3. 

5.3. Temperature and Brightness Profiles 

Figure 4(a) shows the radial profile of ICM temperature, 
based on the result of the spectral fits. The inner 5 annular re- 
gions have the width of about 2.'5, and 20^0% of the detected 
flux comes from the adjacent sky regions. The maximum tem- 
perature within 7' (700 kpc) from the cluster center is ~ 9 keV, 
and it gradually decreases toward the outer region down to ~ 4 
keV around the virial radius. There is a suggestion of weak 
emission with kT ~ 1 keV just outside of r2oo- As shown in 
figure 3(g) and (h), this emission (in magenta) has the intensity 
much weaker than the other components throughout the energy 
range. If one takes into account the systematic error as shown 
in section 5.5, this component does not stay significant in the 
outermost region in figure 3 (h) and in table 5 and 6. Therefore, 
we treat the 1 keV emission in 31.'l-38.'8 to be an upper limit. 

Also, we check projection effect on temperature profiles. 
Based on Fujita et al. (2008), we re-fit each annulus with the 
following method. 

• Fitting the ICM spectrum of the outermost region (24 .'6 - 
3 l.'l ) with a single ICM component. 

• Fitting the next annulus, with a model of a combination 
of the ICM of this annulus and the overlapped emission 
from outer radius, i.e., the best fitting model of the outer- 
most region normalized to account for the spherical pro- 
jection. 

• Proceeding to fitting of inner regions with a model of a 
combination of the ICM of this annulus and the over- 
lapped emission from outer radius, 

The results of this deprojection fitting do not show significant 
difference from the non-deprojection fitting. In sec 6, we use 



the results of non-deprojection fitting method for discussion of 
ICM properties. 

5.4. Electron Density Profile 

The electron density profile was calculated from the normal- 
ization parameter of the Apec model, defined by 

10~ 14 r 

Norm= n e n H dV (1) 

AnD 2 A (l+z) 2 J 

with the unit of crrT 5 , where D A (cm) is the angular diameter 
distance to the source, n e and «h (cm 4 ) are the number densi- 
ties of electron and hydrogen, respectively. We note that the re- 
sultant normalization using an ARF generated by xissimarfgen 
needs the correction by a factor of S OURCE_R ATIO _REG/Q« . 
See Ishisaki et al. 2007 Sec 5.3 for more information. Each 
annular region, projected in the sky, includes emission from 
different densities due to integration along the line of sight. 
We have de-convolved the electron density assuming spheri- 
cal symmetry by successively calculating the value from the 
outermost regions. In this process, we assumed a constant tem- 
perature in each annular region. The uncertainty of the density 
due to this assumption is a few %. 

The resultant radial distribution of n e is shown in figure 4 (b), 
along with the /?-model profile from ROSAT (B = 0.85, r c = 
4.'5). We obtained an upper limit of the electron density in 
the outermost region (3 1 .' 1 — 38 .'8) to be < 2.0 x 10~ 5 crrT 3 . One 
can see that the /3-model gives a fairly good approximation for 
the density profile, even though the temperature shows a large 
deviation from the isothermal case. 

5.5. Systematic Errors 

We have examined the effect of systematic errors on the de- 
rived spectral parameters. We considered 3 components for the 
systematic errors; namely, the NXB intensity with an error of 
+3% (Tawa et al. 2008), fluctuation of the CXB intensity, and 
contamination on the blocking filter. The CXB intensity was 
estimated in section 3.5 with the fluctuation shown in table 3. 
We repeated all the spectral fits by fixing the CXB intensity at 



No. ] 
(a) 0'-2.'5 



X-Ray Study of the Outer Region of Abell 2142 with Suzaku 

(b)2.'5-5.'0 (c)5.'0-7.'8 
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Energy (keV) Energy (keV) Energy (keV) 



Fig. 3. NXB subtracted spectrum for each region. The XIS BI (Black) and FI (Red) spectra are fitted with ICM (wabs + apec), added to CXB+the 
galactic components (LHB, MWH) (apec+wabs(apec+powerlow)). The ICM components are shown in a magenta. The CXB components are shown in 
a black line, and the LHB and MWH emissions are indicated by green and blue lines, respectively. 



the upper and lower boundary values, and the resultant param- 
eters are listed in table 5. The top group shows the parameters 
with the nominal CXB intensity and contamination thickness. 
The parameters in CXBMAX group are for the higher CXB in- 
tensity, and those in CXBMIN are for the lower CXB intensity. 
The effects on the temperature is about 20% in the regions near 
the virial radius. 

The other source of the systematic error is the uncertainty 
in the amount of contamination on the blocking filter of the 
XIS instrument. It showed a gradual increase during the ob- 
servations of A2142, and the standard way is to include 10% 
error in the thickness of the contamination layer. The bottom 
2 groups in tables 5 and 6 show the results. Since this error 
changes the detector response, all the regions are affected by 
the same amount. The change of temperature is 10-20%. 

The Galactic foreground emission is direction dependent 
even though the spectrum can be approximated by the sum of 
two thermal models with temperatures of about 0. 1 keV and 0.3 
keV. We looked into the results for two different backgrounds, 



Offset3 and TCrB regions, which can be treated as the system- 
atic error due to the spatial variation of the Galactic emission. 
As seen in table 6, the effect becomes larger as the position 
goes farther from the cluster center. Near the virial radius, the 
difference in the ICM temperature is about 0.5 keV yielding 
twice different temperatures. 

5. 6. Search for WHIM Lines 

We searched for redshifted O lines in the spectra of all the 
annular regions, by adding gaussian lines in the spectral fit as 
shown in figure 5. The energies of the O VII and O VIII lines, as- 
sumed to have the same redshift of A2142 (0.0909), were fixed 
at 521 eV and 598 eV, respectively. We also assumed the O 
abundance in the ICM to be 0, which gives the highest or most 
conservative upper limits for the O lines. We obtained 2 cr up- 
per limits of the O line intensity for all the regions as shown 
in table 7. The average baryon density in the local universe is 
1.77 x 10- 7 (1 +zf cm" 3 (Takei et al. 2008), or 2.3 x 10~ 7 cirT 3 
at z — 0.0909. Based on the observed upper limits of the line in- 
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Table 5. Best-fit parameters of the ICM with Offset3 background 



Table 6. Best-fit parameters of the ICM with TCrB background 



Rcgioi 



Z(Z ) 



0-2.5 
2.'5-5'0 
5.'0-7.'8 
7^8-10^3 
10'3-12.'9 
r2.'9-16'8 
16.'8-24.'6 
24?6-31.'l 
31.'l-38'8 



-2.5 
5-5'0 

0- 7.'8 

5- ]0'3 
3-12^9 
9-16'8 
8-24^6 

6- 31.'] 

1- 38^8 



0'-2.'5 
2.'5-5.'0 
5.'0-7.'8 
7fS-]0f3 
10'3-12.'9 
12.'9-16'S 
16.'8-24.'6 
24'6-31.'l 
31.'l-38'8 



2.' 5 
5-5'0 
0-7'8 

5- 10'3 
3-12'9 
9-16'8 
8 -24.' 6 

6- 31.'l 



8 

-DIP 

,02+O.W 

9 50+^ 3 ^ 
v - -0.38 
8 47+0-65 

-8-45 

7.79; 

6.14 

3.42+ 

1.01 

1 08+ 5: ^ 
-0-37 



- + m 
in 



0.33+] 
0.28+] 
0.24+j ( 
0.3 1+£ 

019 -0.12 
0.2 (fix) 

0.2 (fix) 

0.2 (fix) 

0.2 (fix) 



F 

^04-IO^V 



A- Z /d.o.f 



L82±0.01 X 10 3 
3.90 + 0.01 X 10 4 
7.50 + 0.01 X 10 3 
2.09 + 0.07 X10 3 
8.30 + 0.32X 10 2 
6.70 + 0.40X 10 2 
1.35 + 0.22X 10 2 
2.00 + 0.11 X10 1 
1.00 + 0.50 



7.74 ± 0.30 X 10 4 
1.65±0.05xl0 4 
3.33±0.O5xl0 3 
9.59 + 0.33 X 10 2 
3.90 + 0.33 X 10 2 
2.67 + 0.22X 10 2 
5.18 + 0.22X 10 1 
6.40 + 0.20 
2.75 + 0.20 



CXBMAX 



TTD-TTT- 



816 L, 

8.99+] 
9.39+] 
8.36+] 
7.65+] 
6.26+ 
3.96+ : _ 
1.06+0f 
1 .02" 



"-17 



33+0.U2 
"■-"-0.02 
28+0.02 

02 3+8:8? 

0.30+|| 

0.2 (fix) 
0.2 (fix) 
0.2 (fix) 
0.2 (fix) 



1.76±0.01xl0 : ' 
3.80±0.05xl0 4 
7.37+0.13x 10 3 
2.05 ± 0.06 X 10 3 
8.18 + 0.32X 10 2 
6.96 + 0.48X 10 2 
1.45 + 0.24X 10 2 
2.16+1.21 X 10 1 
6.80 + 5.80 



CXBMIN 



0'-2.'5 


8 17+0-10 


0.33+0.02 


2.'5-5.'0 


907 +8:l§ 

'■"-0.19 


n 7S +8:81 

0.28_ 002 


5.'0-7.'8 


9 65+ - 38 

-0-38 

8 82 +0.70 


024+0.05 


7^8-10^3 


o 31 +8:8§ 


10'3-12.'9 


8.34+ o: ^ 

629 +?:?5 


19+ o: ^ 
-0.13 


12f9-16'8 


0.2 (fix) 


16.'8-24.'6 


3 3 8 +8:ii 

"lm 


0.2 (fix) 


24'6-31.'l 
31'l-38^8 


1.08+8| 


0.2 (fix) 
0.2 (fix) 



1.87±0.01 X lO 3 
4.04 + 0.05 X10 4 
7.75 + 0. 13x 10 3 
2.19 + 0.06X 10 3 
8.82 + 0.30X 10 2 
7.69 + 0.48 X10 2 
3.47 + 0.25 X 10 2 
1.97 + 1. 20X 10 1 
7.20 + 5.20 



CONTAMI add 



7.85+ 
8.57+ 
9.03+ 
8.18+ 
7.57+j 

5 - 72 -i. 
3.22+ 1 - 



+0.02 

3 



0.33_ 

0.28+] 

0.23 +' 

0.30+^ 

n ih+^ 
0.18_ Q12 

0.2 (fix) 

0.2 (fix) 

0.2 (fix) 



1.86 + 0.01 x 10 3 
3.94 + 0.05 xlO 4 
7.62 + 0. 13x 10 3 
2.15 + 0.06X 10 3 
8.50 + 0.03 XlO 2 
6.96 + 0.48X 10 2 
1.45 + 0.24X 10 1 
l.OSiO^OxlO 1 



0.92 
0.93 
1.03 
1.07 
1.19 
0.85 
0.92 
0.94 



CONTAMI Wk red 



9.22+ - 19 

981 +o:38 

8.80+H 
8.39+; 
6.60+!;" 
415 +3.52 



+0.02 



0.34; 

0.28+] 

0.24+] 

0.31+] 

0.18+] 
0.2 (fix) 
0.2 (fix) 
0.2 (fix) 
0.2 (fix) 
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1.81 ±0.01 x 10^ 
3.89 ± 0.05 X 10 4 
7.50 + 0.13X 10 3 
2.09 + 0.06 X 10 3 
8.18 + 0.32X 10 2 
6.56 + 0.40X 10 2 
1.24 + 0.24X 10 2 
1.91 + 1. 10X10 1 
9.10 + 5.40 



7.71+0.03X10 4 
1.64 + 0.05 XlO 4 
3.31 +0.06X 10 3 
9.54 + 0.34X 10 2 
3.86 + 0.34X 10 2 
2.60 + 0.20X 10 2 
4.89 + 0.18X 10 1 
7.00 ±0.20 
3.15 + 0.03 



0.93 
0.94 
1.04 
1.07 
1.19 
0.85 
0.93 
0.94 
0.87 



712/604 
691 /604 
681 /604 
309/270 
338/270 
381/358 
272/217 
238 / 223 
299 / 223 



7.52 + 0.30X 10 4 


0.93 


712/604 


1.60 + 0.06X 10 4 


0.93 


690 / 604 


3.23 + 0.O5X10 3 


1.03 


675 / 604 


9.35 + 0.32X 10 2 


1.07 


301 /270 


3.82 + 0.32X 10 2 


1.18 


337/270 


2.75 + 0.23 X 10 2 


0.84 


376/358 


5.64 ± 0.29 Xl0 ] 


0.90 


273/217 


7.55 + 0.40 


0.89 


231 /223 


2.35 + 0.03 


0.82 


288 / 223 



7.98 + 0.31 x 10 4 


0.92 


713/604 


1.70 + 0.06X 10 4 


0.93 


693 / 604 


3.44 + 0.05 XlO 3 


1.03 


687 / 604 


1.00 + 0.03 XlO 3 


1.07 


3 1 1 / 270 


4.14 + 0.34X 10 2 


1.18 


339/270 


2.33 ± 1.81 X 10 2 


1.19 


431/358 


7.35 + 0.60X 10 1 


1.26 


255/217 


7.05 + 0.20 


0.97 


242 / 223 


2.55 ±0.02 


0.94 


304 / 223 



7.86±0.03xl0 4 


0.92 


734 / 604 


1.67±0.06x 10 4 


0.93 


710/604 


3.73±0.O5xl0 3 


1.03 


685/604 


9.74±0.29x 10 2 


1.06 


300/270 


3.97±0.30x 10 2 


1.18 


325/270 


2.73±0.23xl0 2 


0.84 


376/358 


5.48±0.27x 10 1 


0.91 


264/217 


4.00 ±1.00x10° 


0.95 


255/223 




1:0.88 


319/223 



764 / 604 
747/604 
706/604 
322 / 270 
357/270 
392/358 
286/217 
226 / 223 
283/223 



*: Normalization of the apec component scaled with a factor SOURCE_RATIO_REG f>e from lablo , 
N o rm= SOURCE.RATIOJ!EG fnen H dV/(4»(l+z^)D^) X 10-^cm-5 arcmh,- 2 , 
where is the angular diameter distance to the source, 
t: 10" 8 phoionscm" 2 s _1 arcmin" 2 . Energy band is 0.4- 10.0 keV. 

Surface brightness of the apec component scaled with a factor SOURCE-RATIO-REGjQt from lable 3. 



Region 


kT(keV) 


Z(Zo) 


norm* 


''0 4-lO^V 


BI/FI 


* 2 /d.o.f 


0'-2.'5 


8 16+010 


033+0.02 
°- 28 -0.02 


1.82 + 0.01 X 10 3 


7.74±0.03x 10 4 


0.92 


711/604 


2.'5-5'0 


O02 +8:l§ 

-0.19 


3.89 ±0.05 XlO 4 


1.65±0.06x 10 4 


0.93 


689/604 


5.'0-7?8 


n 47+0.38 
* -0.38 


02 4+0.05 


7.50±0.13xl0 3 


3.33 ±0.04 XlO 3 


1.03 


674/604 


7.'8-10.'3 


847 +0.66 


0.31+^ 


2.12±O.07xlO 3 


9.63 ±0.29 XlO 2 


1.06 


298 / 270 


10.'3-12'9 


784 +8:M 


19+^ ^ 
UIV -0.12 


8.39 ±0.32 XlO 2 


3.94 ±0.31 XlO 2 


1.17 


323/270 


12.'9-16'8 


6.02+^ 


0.2 (fix) 


6.88 ±0.48 XlO 2 


2.70 ±0.24 XlO 2 


0.84 


373/358 


16'8-24.'6 
24.'6-3l'l 


3.6 2 +?:7 l 6 
107 +o:5! 

1 07+ o; ^ 
'■"'-0 30 


0.2 (fix) 
0.2 (fix) 


1.43 ±0.24 XlO 2 
1.91 ± 1.21 X 10 1 


5.47 ± 3.05 X 10 1 
6.75 ±0.25 


0.89 
0.97 


264/217 
228/223 


3lf 1-38.'S 


0.2 (fix) 


9.11 ±6.07 


3.15 + 0.25 


0.87 


275 / 223 


CXBMAX 


0'-2.'5 


g .l 6 +lUO 


33+0.02 

-0.02 

28+0 ° 2 
024+8$ 


1.82 + 0.01 x 10 3 


7.73 ±0.03 xlO 4 


0.92 


711/604 


2.'5-5'0 


„,+8:l§ 


3.89 ±0.05 XlO 4 


1.65±0.06x 10 4 


0.93 


689/604 


5.'0-7?8 


946+8:3! 
g 46 +8:S 


7.50±0.13xl0 3 


3.31±O.05xlO 3 


1.03 


674/604 


7.'8-10.'3 


0.31+^ 


2.09 ±0.76 XlO 3 


9.54 ±0.30 XlO 2 


1.06 


299 / 270 


10.'3-12'9 




19+^ ^ 
-0.13 


8.18±0.32xl0 2 


3.86±0.31xl0 2 


1.18 


324 / 270 


12,'9-16'8 




0.2 (fix) 


6.56 ±0.40 XlO 2 


2.57 ±0.22 XlO 2 


0.84 


373 / 358 


16.'8-24'6 


":;+?:§§ 


0.2 (fix) 


1.28±0.24xl0 1 


4.81 ±0.22x10* 


0.91 


267/217 


24.'6-3l'l 


2 10+° : ^ 
Z - 1U -0.32 


0.2 (fix) 


7.00 ±1.55 


2.55 + 0.15 


0.91 


236/223 


3lM38.'8 










0.87 


287/223 



0'-2.'5 
2.' 5 -5^0 
5.' 0-7^ 8 
7.'8-10.'3 
10^3-12.'9 
12.9-16.'8 
16' 8 -24.' 6 
24 ^ 6-31.' 1 
3l'l-38.'8 



0'-2.'5 


8.16 


2.' 5 -5^0 


9.03 


5.' 0-7^ 8 


9.51 


7.' 8- 10.' 3 


8.53 


10^3-12.'9 


7.95 


12.9-16.'8 


6.57 


16.' 8 -24.' 6 


3.96 


24. 6-31.' 1 


1.15 


3l'l-38.'8 


1.09 



+0.10 
-0.10 
+0.19 

+8:4? 



+?:?? 

-1.12 
,+2.28 

M 



33+11.02 

UJJ -0.02 


1.82 + 0.01 x 10 3 


7.74±0.30x 10 4 


0.92 


712/604 


28+0.02 

024+8:8? 

-0.05 


3.89 + 0.05 XlO 4 


1.65 ±0.06 x 10 4 


0.93 


689/604 


7.62 + 0.13X10 3 


3.34 + 0.51 x 10 3 


1.03 


675 / 604 


„.3 ,+0.09 


2.12 + 0.07X10 3 


9.69 + 0.36 XlO 2 


1.07 


299/270 


, 9 +8:?? 

-0.12 


8.50 + 0.32X10 2 


3.98 + 0.32 XlO 2 


1.18 


324/270 


0.2 (fix) 


7.28 + 0.48X10 2 


2.79 + 0.24 XlO 2 


0.84 


372/358 


0.2 (fix) 


1.47 + 0.22 XlO 2 


5.71+0.28X10 1 


0.91 


268/217 


0.2 (Ox) 


2.23 + 0.14X 10 1 


8.10 + 0.60 


0.99 


222/223 


0.2 (fix) 


9.57 + 6.30 


3.35 + 0.15 


0.92 


273 / 223 



CONTAMI 10',-i ADD 



0'-2.'5 


7 86 +o.io 


0.33 +(J -^ 
28+^ 


1.86±0.01 X 10 3 


7.86 ± 0.32 X 10 4 


0.92 


734/604 


2.' 5 -5^0 


8 57+8:28 


3.94±0.05xl0 4 


1.67 ±0.06 X 10 4 


0.93 


710/604 


5.' 0-7^ 8 


9 05+°^ 
8.21_ 047 


023 +8:Bi 


7.62±0.13xl0 3 


3.37 ±0.05 XlO 3 


1.03 


684 / 604 


7.'8-10.'3 


30+^ 
U " 5U -0.O8 


2.15±0.66xl0 3 


9.73 ±0.25 XlO 2 


1.06 


299 / 270 


10^3-12.'9 
12^9-16.'8 


764 +0.80 

5 83+i:I8 


019 +0.12 
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^Normalization of the apec component scaled with a factor SOURCE_RATIO_REG Q e from lable , 
Norm= SOURCE_RATIO.REG j^dv/^g + 2 2 )p 2 , x 1Q -22 cm -5 arcmin -2 _ 

where Z)y, is llie angular diaineler dlslanee lo ihe source. 

t: 1 - ^ pholonscm -2 s~ ^ arcmin -2 . Energy band is 0.4 - 10.0 keV. 

Surface brightness of Ihe apec component scaled with a factor SOURCE-RATIO-REGi2 e from table 3. 



tensities, we derived upper limits of the over-density (6) of the 
gas assuming the line-of-sight depth of 2 Mpc and the temper- 
ature of 2 x 10 6 K. These assumptions are the same as in Takei 
et al. (2008). The results are summarized in table 7 for the 
two different Galactic backgrounds separately. We note that the 
lowest upper limits, with overdensity less than 280 (assuming 
background of Offset3) or 380 (assuming background of TCrB) 
is obtained for the outermost region (r = 3 1 '.0 — 38.7 ). We also 
derived upper limits of the overdensity from the measured elec- 
tron density. In sec 5.4, we showed the 1 <r upper-limit of the 
electron density to be <2.0x 10~ 5 crrT 3 in the outermost region. 
In the range of 2 <x, the upper limit of the electron density is 
< 2.6 x 10" 5 cm 4 , corresponding to an overdensity < 113. 



6. Discussion 

Suzaku performed four pointing observations in A2142 and 
its outside regions along the merger axis. The temperature was 
measured out to the virial radius (r2oo ~ 2.5 Mpc) for the first 
time. The ICM temperature was found to drop from about 9 
keV around the center to about 3 keV at r2oo- We detected no 
significant signal from the WHIM and set upper limits of its 
over-density. We will attempt to evaluate the cluster properties 
(temperature, electron density, and entropy) and discuss their 
implications. 

6.1. Temperature Profile 

Some numerical simulations (Ettori et al. 2004; Borgani et 
al. 2004) predicted that the intracluster gas temperature drops 
to about 50% of the central temperature around r2oo under 
hydrostatic equilibrium. Those results reproduced observed 
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Fig. 4. Radial profiles for (a) temperature, (b) 3-dimensional electron 
density. In the temperature plot, the uncertainty range due to the combi- 
nation of ± 3% variation of the NXB level and the maximum/minimum 
fluctuation of CXB is shown by two green dashed lines. In a similar 
way, the uncertainty range due to the effect of contamination on the 
blocking filter is shown by red lines. Orange diamonds show Chandra 
data from Markevitch et al. (2000). 

temperature profiles to about 0.5r2oo (De Grandi & Molendi 
2002; Vikhlinin et al. 2005; Pratt et al. 2007). In particular, 
recent XMM-Newton (Pratt et al. 2007) results showed tem- 
perature profiles out to 0.8r2oo- 

We compare here the observed temperature profile of A2142 
with the pro?le of other clusters that have been observed with 
Suzaku. Table 8 shows the list of clusters for which Suzaku 
measured the temperature out to r2oo- Their temperature pro- 
files are shown in figure 6. All the temperatures are normal- 
ized by the flux-weighted average of each clusters. The results 
clearly show a systematic drop of the temperature by a factor 
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Fig. 5. Fits for the constraint on the intensities of the /ovn and /ovm 
lines in the 31.'l-38'8 region. The /ovn and /ovm emission lines are 
shown as a cyan lines, and the notations of the other lines are same as 
shown in figure 3. 

of 3-5 from the center to r2oo- 

Among these clusters, A1689 shows a significant directional 
difference in the sense that the temperature drop is small along 
the filament direction (Kawaharada et al. 2010). Another sys- 
tem A1413 also indicates a somewhat flatter temperature pro- 
file, and the measured direction is along the longer axis of its 
X-ray elongation (Hoshino et al. 2010). 

Burns et al. (2010) discuss a non-equilibrium effect on the 
temperature profile based on N-body + hydrodynamic simula- 
tions. Those simulations indicate cluster temperature decline 
by factor of ~ 3 at r2oo, consistent with the feature observed for 
two nearby clusters, PKS 0745-191 and A 1795 (George et al. 
2008; Bautz et al. 2009). Burns et al. (2010) approximated the 
average temperature profile by a function, 



= A 



l+B 



(2) 



They obtained the best-fit values as A = 1.74 + 0.03, B = 
0.64 + 0.10, B = -3.2 + 0.4 for the two cluster data. In fig- 
ure 6, we plotted model prediction temperature profile. Black 
dotted curve shows the best-fit temperature profile, and dashed 
lines show lcr error range reported by Burns et al. (2010). The 
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relation of Burns et al. (2010) represents the temperature pro- 
files for the 6 clusters fairly well. This approximate "universal" 
temperature profile suggests that clusters generally hold self- 
similar relation even near r2oo, where some temporary effects 
caused by infalling matter may be seen. Note that the tem- 
perature drop in the filament direction of A 1689 (not shown 
in Fig 6) is flatter than the average A 1689 profile, suggesting 
a very efficient heating going on in the filament direction as 
compared with typical clusters. 

Table 8. Cluster samples and those properties 



Cluster 




Ref. 


k(T) 


loo 








[keV] 


[Mpc](arcmin) 


Abell 2142 


0.090 


This work 


8.6 


2.46 (24.6) 
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George et al. 2008 
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2.21 (19.6) 


Abell 2052 


0.036 


Tamura et al. 2008 


3.2 


1.54 (36.7) 
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Reiprich et al. 2009 
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Fig. 6. Scaled projected temperature profiles. The profiles have been 
normalized to the mean temperature. The r2oo derived from Henry et 
al. (2009). Dotted line show simulation result Burns et al. (2010). Two 
gray dashed lines show standard deviation. 

All the clusters in figure 6 are morphologically relaxed sys- 
tems, even though all of them show an elliptical shape to some 
extent. Therefore, we may regard the model temperature pro- 
file to be a typical feature for clusters having relaxed morpholo- 
gies. We note that the present A2142 profile shows a smooth 
temperature decline with radius, even though it is a merger 
cluster characterized by clear cold fronts within 3' (0.11r2oo) 
from the center. 

We note that when we fit the temperature profiles with 
the above formula 2, A2142 data are characterized by (3 — 
-5.8 + 0.8, which indicates a steepr slpe than the average, 
/3 = -3.2 + 0.4. On the other hand in A1689, the filament di- 
rection shows a flatter temperature profile (Kawaharada et al. 
2010) than the average curve. It has been discussed that A 1689 
still holds the heating feature caused by recent matter infall. 
The same process may be working in A1413 which also shows 



somewhat flatter temperature distribution (Hoshino et al. 2010) 
even though the presence of filament is not clear in the ob- 
served direction. 

Even though flatter temperature profiles are regarded as a 
sign of recent matter infall, steeper temperature curves such 
as seen in A 2142 and abrupt temperature drops, as seen in 
A 3667 and Virgo Cluster can also be a sign of recent heat- 
ing (Finoguenov et al. 2010; Urban et al. 2011). This point 
will be discussed further in the next section. Therefore, simple 
steepness of the temperature profile seems to carry somewhat 
degenerate information about the cluster evolution. 

6.2. Entropy Profiles 

The entropy of ICM is used as an indicator of the energy ac- 
quired by the gas. We will refer to "entropy" of the ICM by 

—2/3 

K = kTn e following the recent convention. Numerical simu- 
lations indicate that a self-similar growth of clusters commonly 
show entropy profiles approximated by r up to r2oo, exclud- 
ing the cool core region (Voit et al. 2003). Recent XMM- 
Newton results on the entropy profiles of 3 1 clusters showed 
a median slope of 0.98 out to about rsoo which is approxi- 
mately 0.5r2oo (Pratt et al. 2010). The slope also showed posi- 
tive correlation with the average temperature. They also found 
that morphologically disturbed clusters show a large scatter 
(0.5-2.0) in the slope. Suzaku has extended the entropy mea- 
surement close to r2oo for several clusters, and showed a flat- 
tening or even a decrease at r > 0.5r2oo as mentioned earlier. 
Discussion has been made that the ICM may not be under equi- 
librium in the cluster outer regions (Hoshino et al. 2010). 

Figure 7 shows the entropy profile of A2142 based on the 
present Suzaku measurements. The solid line indicates the 
slope of 1.1. The entropy slope is consistent with this value 
in the radius range 0.1 — 0.4r2oo (2'5 - 10')- The slope becomes 
flatter at r > 0.4r2oo and goes negative near ^200- This feature is 
more clearly seen in figure 7(b), which shows the ratio to the 
r curve which can be fitted for the inner profile, r < 0.4r2oo- 
Including the previous results for PKS0745-191 and A1413, 
flattening of the entropy profile in such outer regions looks to 
be a common feature. 

As for the cause of the entropy drop relative to the predicted 
r • relation, two possibilities have been pointed out so far. One 
is the non-equilibrium effect (Hoshino et al. 2010; Burns et al. 
2010) and the other is the dumpiness (Simionescu et al. 201 1), 
both working only in the outer regions. These two conditions 
are not mutually exclusive but may be closely related with each 
other. 

Regarding the effect of dumpiness, Nagai (2011) exam- 
ined with numerical simulations how much clumping in ac- 
creting gas could cause overestimation of gas density since the 
emissivity scales as density squared. This density overestima- 
tion gives lower entropy values, hence resulting in a flatten- 
ing around r2oo- Results of Nagai (2011) show that the en- 
tropy drops by 10-15% from the true value at around r2oo- In 
the A2142 case, the observed entropy curve starts to flatten 
around rsoo and decreases to less than 40% of the r l 1 extention 
at ?"2oo (~ 25')- This large suppression of the entropy seems 
too large according to the simulation. Therefore, clumping is 
not likely to be the only or the dominant source of the entropy 
suppression. 
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Fig. 7. (a) Entropy profile: Diamonds show Suzaku results, and solid 
straight line shows the universal trend of the entropy curve ~ r 1,1 , re- 
ported by Voit (2005). (b) Entropy ratio relative to the r u profile. Gray 
diamonds show the result of A1413 (Hoshino et al. 2010) (c) t„ (dia- 
monds) as a function of radius, compared with estimated t aiipscli (solid 
line). Different lines show different shock speeds (gray dash:400 km 
s~\ dash: 800 km s" 1 , dot: 1200 km s" 1 ). 



6.3. Ion-Electron Relaxation 

The entropy profile of A2142 shows a flatter slope than the 
predicted one based on numerical simulations, and even a neg- 
ative slope is seen in the outermost region. Kawaharada et al. 
(2010) showed a flat entropy profile in the filament direction of 
A1689 and discussed effects from surrounding environments. 



If we make a simple assumption that the entropy profile of ICM 
should follow r 1 1 (Voit et al. 2003) under an equilibrium, the 
present result for A2142 suggests that electrons in the cluster 
outer region is not in the equilibrium. 

Hoshino et al. (2010) discussed the possible difference be- 
tween ion and electron temperatures in the outer low-density 
region, where the relaxation time scale is inversely proportional 
to particle density. We follow their estimation and look into 
the behavior of ion (r,) and electron (T e ) temperatures. The 
following assumptions are made. 

1. Ions are initially heated through accretion shocks at r2oo- 

2. Ions achieve thermal equilibrium with a timescale of ion- 
ion relaxation, f„, after this heating. 

3. Thermal energy is transferred from ions to electrons 
through Coulomb collisions. This process takes a 
timescale of t ei which is much longer than both f„ and 
the electron-electron relaxation time t ee . 

The relation between T, and T e can be estimated in the 
following way after Hoshino et al. (2010). We look into a 
position-dependent time, f e i apse d, which is a rough measure of 
the time elapsed from the shock heating. We assume that the 
shock heating takes place at the virial radius and that the shock 
wave propagates with a constant speed (400, 800, 1200 km 
s _1 ) in ICM toward inner regions. This f e i apse d can be compared 
with the equilibration timescale t ei , and we may assume that 
electrons are fully heated up when f e iapsed > t„. In figure 7(c), 
we compare f e i apse d wltn t e i as a function of radius. Each curves 
show the position-dependent f e i apse d, assuming that the gas falls 
through either free-fall or constant velocities. The equilibration 
timescale t ei is significantly longer than f e i apse d in the cluster 
outskirts at r > 20', suggesting that the ion temperature can be 
higher than the electron temperature. This feature is essentially 
the same in A 1413 as reported by Hoshino et al. (2010). 

Theoretical studies of non-equilibrium ionization state and 
an electron-ion two-temperature structure of ICM in merg- 
ing galaxy clusters have been carried out (Rudd & Nagai 
2009; Akahori & Yoshikawa 2010). They show that in merging 
clusters, T e is lower than the average temperature by 20-30%. 
Following Hoshino et al. (2010), we estimate the possible de- 
viation of T e by assuming that the average temperature is given 
by the entropy and density as fcr gas = Kn 2 ^, and further that 
the entropy K follows the power-law profile oc r 1 1 to r2oo (Voit 
et al. 2003). Figure 7(b) shows thus estimated ratio r c /r gas . 
Previous A1413 results by Hoshino et al. (2010) are shown 
for compariton. A2142 values are lower than A1413 at r2oo, 
which may reflect stronger suppresion of T e in merging clus- 
ters. The r e /r gas value in A2142 is 0.51+j}^ at r 200 , which 
is in good agreement with the theoretical result for unrelaxed 
clusters (Rudd & Nagai 2009). 

Even though this feature suggests that T e is substantially 
lower than the average gas temperature around r2oo, we have 
to note that the temperature profile of A2142 agrees with those 
for other relaxed clusters as seen in figure 6. This may sug- 
gest that the entropy drop is rather due to high n e around r2oo- 
However, as shown in figure 4 (b), the observed n e profile also 
follows a smooth /3-model and no peculiar hump is seen in the 
outer region. It may be that the outer region of A2142 has 
not experienced strong merger recently, but then the entropy 
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Our differential mass density profile is consistent with the pre- 
vious weak-lensing result which shows c = 4.26*^ (Umetsu 
etal. 2009). 

We can evaluate the contribution of non-thermal pressure by 
comparing the X-ray and the weak-lensing masses. The X- 
ray to lensing mass ratio is 0.90^q jj, which indicates that, by 
taking the lowest boundary value, the maximum non-thermal 
pressure can be 45%. As shown in §6.3, we discussed the pos- 
sible difference between ion and electron temperatures as the 
cause of the entropy drop from the r u relation. In order for the 
entropy profile to follow the r l i relation after electron temper- 
ature reaches the ion temperature, ion temperature should have 
higher temperature than electrons. The ratio of electron to ion 
temperature is estimated to be 0.5 ^ around the virial radius 
(see fig7b: this value also indicates the ratio of electron to ion 
temperatures). Then, since the gravitational force is balanced 
with the sum of electron and ion pressure in hydrostatic equi- 
librium, the gravitational mass would rise roughly by the same 
factor. Although this leads to poorer match between X-ray and 
lensing mass, they are still barely consistent within statistical 
and systematic uncertainties. 

We note that non-thermal pressure can make additional con- 
tribution to the cluster mass estimation. Dynamical effects such 
as turbulence and bulk motion of ICM are not yet measured, 
but can give significant effects in the cluster outer regions. 
Some numerical simulations (Nagai et al. 2007; Piffaretti & 
Valdarnini 2008) predicted that such non-thermal pressure 
could add up 15-30% of the cluster mass. 

We also looked into the gas mass distribution using the ob- 
tained electron density profile and showed the results in fig- 
ure 8(a) with black crosses. Here, the gas mass does not in- 
clude the stellar mass. The resultant gas mass fraction at r200 
is \A.A + J 4 \%. This agrees with the expected hot-gas fraction 
(15%) in the universe (Komatsu et al. 2011), and consistent 
with the previous result, ~ 18% at 22'. 8, incorporating the weak 
lensing and Sunyaev and Zel'dovich observations (Umetsu et 
al. 2009). This indicates that, in A2142, there is no strong need 
to invoke the gas dumpiness to account for the gas mass frac- 
tion as in the Perseus cluster (Simionescu et al. 201 1). 



"saturation" in the outer region will have to have a more uni- 
versal origin. We certainly need to look into this problem with 
a wider range of cluster sample regarding the cluster size and 
morphology. 

6.4. Electron Density Profile 

The electron number density profile in figure 4 shows a de- 
crease down to ~ 10~ 5 cirT 3 around the virial radius. ROSAT 
study showed that the electron density profile could be fitted 
well except for the cluster center with the /?-model, after mod- 
ification for the cluster ellipticity (Henry & Briel 1996). In the 
outskirts of relaxed clusters, many systems indicate a density 
profile of n e oc r 22 (Zhang et al. 2006). Also, n e oc r -i-8±o.28 is 
obtained for the REXCESS sample (Croston et al. 2008) 

The measured A2142 density profile is fitted with a power- 
law model with an index of -1.98 + 0.13 for the entire radial 
range and -2.53+0.25 in the outskirts only (r > 0.5r 20 o). These 
results agree with the previous studies for other clusters, and 
with the predicted profile of n e oc r~ 2 - 5 based on the yS-model 
(J3 = 0.85) in the outer region. 

Recently, Eckert et al. (2011) pointed out a factor of about 
3 discrepancy in the surface brightness in the outer region of 
PKS 0745-191 between Suzaku and ROSAT results. They ar- 
gue that it is most likely to be caused by an incorrect sub- 
traction of the Galactic emission for the Suzaku data. As for 
A2142, ROSAT PSPC observed its emission to r 20 o (Henry & 
Briel 1996) and showed n e ~5x 10~ 5 cm~ 3 at r2oo- The present 
n e profile, shown in figure 4, indicates a lower but consistent 
value with the ROSAT one within a factor of ~ 1.5. Therefore, 
our subtraction of the Galactic background in the nearby sky 
region is considered to give reliable results. 

6.5. Mass Estimation to r2oo 

We estimate the gravitational mass of A2142 to r2oo based 
on the observed temperature and density profiles. Here, we 
assume hydrostatic equilibrium and spherical symmetry, and 
calculate the gravitational mass within 3-demensional radius R 
with the following formula (Fabricant et al. 1980), 

R 2 dP g kTR Idlnpg d\nT\ 

Mr- = + , (3) 

PgG dR p-m p G \ dlnR d\\\R j 

where G is the gravitational constant, and yu(« 0.6) is the mean 
molecular weight of gas. 

In figure 8, solid diamonds show the gravitational mass of 
A2142 based on the observed temperature and density pro- 
files. Two solid vertical lines (black and gray) indicate the 
mass derived from the previous weak lensing studies which 
are summarized in table 9 (Okabe & Umetsu 2008; Umetsu 
et al. 2009). Dashed-line diamonds show gravitational mass 
using the y6-model density profile together with the observed 
temperature gradient. The resultant mass agrees well with the 
previous studies around the virial radius. 

We compare the parameters of an NFW profiles, defined as 



(r/r s )(l+r/r s ) 2 

where p s is central density parameter and r s is the scaled ra- 
dius. Figure 8(b) show the resultant values of differential mass 
density (c = 2.8^ 2), which reflect the dark matter potential. 



Table 9. Mass estimation of A2142 

Reference r M200 

arcmin 10 14 M o 

Okabe & Umetsu 2008 722 13.7 + 6.0 

Umetsu et al. 2009 22.2 12.3!?,;° 

This work (Observed gas density profile) 22.5+3.9 11. I*? - ? 

This work (yS-model gas density profile) 22.5+3.9 10.4^ 



7. Summary 

We observed Abell 2142 in the direction of the possible 
merger axis with Suzaku and detected the ICM emission up to 
the virial radius r2oo (2.5 Mpc ~ 24 .'6). We derived radial pro- 
files of temperature, electron density, and entropy, and com- 
pared these properties with the previous results for other re- 
laxed clusters. We summarize the main features of A2142 as 
follows; 
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Fig. 8. (a) Total integrated gravitational mass M<« profile of A2142. Black diamonds show gravitational mass profile estimated from the Suzaku data. 
Gray dashed diamonds show a profile assuming the electron density to follow fi model distribution with /? = 0.85. Black and gray vertical lines show 
masses estimated by weak-lensing analysis summarized in table 9 (Okabe & Umetsu 2008; Umetsu et al. 2009). Black cross show gas mass profile 
estimated from the observed electron density profiles, (b) Same as (a), but for differential mass density. We ignore the data around 10', because it gives 
an unphysical negative value, Dotted curve show NFW-model results by Umetsu et al. (2009). Two dashed lines show standard deviation. 



• The ICM temperature gradually decreases toward the 
outer region from about 10 keV at 0.2r2oo to about 4 keV 
at r 2Q Q. 

• The temperature profile in the outer region of A2142 
agrees well with the results of other clusters observed 
by Suzaku. 

• The average temperature profile for different clusters can 
be described by the formula by Burns et al. (2010) up 
to r2oo, suggesting that non-thermal pressure support is 
significant in the outer regions. 

• The electron density profile decreases down to ~ 
10~ 5 cirT 3 at the virial radius and well agrees with the 
£ model with /? = 0.85. 

• The entropy profile within about 0.4r2oo follows r 11 , pre- 
dicted by the accretion shock heating model. The profile 
becomes flatter and finally shows a negative slope around 
r2oo, suggesting significant deviation from the equilib- 
rium condition. 

• Based on the temperature and entropy profiles, and re- 
quired relaxation time, we discuss that T, is likely higher 
than T e in the outer region of A2142. 

• The derived mass profile is in agreement with the weak- 
lensing mass. The difference between X-ray and lensing 
masses allowed by the error can explain the contribution 
of the non-equilibrium effect. 

Resent Suzaku results on the temperature and entropy mea- 
surements to r2oo, including the A2142 ones, jointly suggest 
that ICM in cluster outskirts is deviating from thermal equilib- 
rium in the sense that the electron temperature could be signif- 
icantly lower than the ion temperature. This suggests that the 
ICM around r2oo has experienced bulk motions and/or turbu- 
lence within a time scale of about 10 9 yr. Such gas motions 



with a velocity of a few hundred to thousand km s _1 can be ob- 
servable in future by X-ray microcalorimeters, such as by SXS 
instrument on ASTRO-H (Mitsuda et al. 2010) which has 20- 
30 times higher energy resolution than CCD instruments. Since 
the FOV of SXS is rather small (3' x 3'), expected counts from 
Fe-K lines are typically around 10 and 60 at r2oo and 0.5r2oo, 
respectively, for a 300 ksec exposure. Therefore, long expo- 
sures for selected good targets with ASTRO-H will be needed 
to show the dynamical features appearing in the cluster outer 
regions for the first time. 

We are grateful to N. Okabe and K. Umetsu for informa- 
tion about thier Suzaku data and N. Okabe providing SDSS 
galaxy distribution data. H. A. was supported by Grant-in- Aid 
for JSPS Fellows (22-1582) and the MEXT program " Support 
Program for Improving Graduate School Education". NO ac- 
knowledges financial support from Grant-in- Aid for Scientific 
Research No.22740124. 

Appendix 1. Point Source Analysis 

As for the point-source subtraction, we first ana- 
lyzed the XMM-Newton data (Observation ID=01 11870101, 
01 1 1870401) in which faint sources were resolved better than 
the Suzaku data. The data covered up to the virial radius (~ 2.5 
Mpc). We used wavdetect tool in CIAO (CIAO version:4.0.1) 
to detect point sources and used surround annular region for 
background subtraction. We summed MOS 1 and MOS2 spec- 
tra to increase statistics, and fitted by pegpowerlaw model 
which offered photon index and flux in selected energy band. 
In the outer region (r > 25.'6), we removed the sources from 
the Suzaku data which were selected by eye. In the XMM- 
Newton case, the source extraction radius is 30", and the sur- 
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Table 10. Informations of point source in two XMM and Suzaku observations of A2142 center, OFFSET1, OFFSET2 OFFSET3. 



[Vol. 



XMM-Newton(MOSl+MOS2) Suzaku 



No.* 


Photon Index 


Flux* 


* 2 M.o.f 


Photon Index 


Flux* 


^ 2 /d.o.f 


1 (239.401,27.485) 


1 




108.9/63 


1 96 +(uy 

-0.23 


1 02 +uj; 

-Q-3Q 


53.0/52 


2 (239.393, 27.287) 


31 +0 5f ' 

u - ,i -0.31 


35.6/23 




1 40 +ilS 

1 - HU -0.48 


59.0/53 


3 (239.334, 27.275) 


9 9 1 +8:tSf 

Z ' -0.50 


i 97+0.46 
1,z,/ -0.18 


34.0/31 


2.08!;^ 


1 8 C +0.46 
1.88_q.42 


75.0/53 


4 (239.283, 27.366) 








1 64 +0 - 23 


9 09+0.59 
Z " J -0.59 


74.0/52 


5 (239.532,27.351) 








94+ - 54 

-Q.74 


2 ' 28 -179 


58.0/51 


6 (239.295, 27.605) 








L89+ li.i5 


89 + " : ^ 
-Q4S 

3 ' 23+ 039 


70.0/54 


8 (238.922, 27.686) 








1 qa+oMS 
-Q-1Q 
1 63 +015 


60.0/54 


9 (239.012,27.773) 








1.65^1 


52.0/54 


*:Serial number for point source. 




f :Position of the point 


source. 







t:The 2.0-10.0 keV flux in units of 10~ 13 erg cirT 2 s _1 



rounding background ring in estimating the flux is defined by 
30" -60", respectively. In the Suzaku case, those are 1' and 
1' - 2' , respectively. We found 3 sources in the XMM-Newton 
data whose fluxes were higher than 3 x 10~ 14 erg cirT 2 s _1 in 
the energy range 2-10 keV. In the Suzaku case, we found 5 
sources whose fluxes were higher than 8 X 10" 14 erg cirT 2 s _1 
in the energy range 2-10 keV. 
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